Population inversion of driven two-level systems in a structureless bath 
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We derive a master equation for a driven double-dot damped by an unstructured phonon bath, 
and calculate the spectral density. We find that bath mediated photon absorption is important at 
relatively strong driving, and may even dominate the dynamics, inducing population inversion of 
the double dot system. This phenomenon is consistent with recent experimental observations. 
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Quantum dots arc of great current interest due to their 
controlability, their sensitivity to the quantum nature of 
the phonon environment [1, 2, 3], and their potential ap- 
pUcation to quantum information processing [4] . Double- 
dot systems, formed by depleting a 2D electron gas down 
to few conduction band electrons in two wells, have shown 
coherent quantum phenomena [5, 6]. 

For the last decade, sensitive electrometry has been 
able to probe double-dot systems in the single electron 
regime [7, 8]. Microwave spectroscopy reveals peaks 
when the double-dot is resonant with a driving field, in- 
dicating Rabi oscillation [6, 9]. Similar phenomena have 
been seen in driven superconducting systems [10, 11, 12]. 

In this paper we study the interaction of a phonon bath 
with a strongly driven double-dot. Other work on driven, 
dissipative two- level systems (2LSs) is reviewed in [13]. 
We begin by deriving a spin-boson model from micro- 
scopic principles, from which we develop a master equa- 
tion for the interaction between the double-dot and an 
unstructured phonon bath (i.e. one lacking resonances). 
We analyse the model and find that at zero temperature 
and non-resonant driving, phonon assisted photon ab- 
sorption becomes significant. The rate balance between 
this and the ever-present phonon induced relaxation re- 
sults in an increased occupation of the excited state. 

We show that for strong cw-driving and sub- quadratic 
spectral densities, there are regimes where bath-mediated 
excitation dominates, leading to the surprising prediction 
of a large steady-state population inversion in a 2LS. This 
is consistent with recent experimental results [6, 14]. 

While inversion is integral to lasers, three levels are 
usually required. Inversion of nominal 2LSs has been 
discussed for superconducting systems [15, 16], Er-doped 
glass fibre [17] and strongly driven atom-optical systems 
[18, 19, 20]. These focus on structured spectral densities 
or extremely strong driving in contrast with the mecha- 
nism presented here, which only requires significant cou- 
pling to the unstructured bath near the Rabi frequency. 

Model. We consider a periodically driven 2LS (a 
qubit), modulating both bias and tunnelling, 

Hq = —{t<Jz + /S.ax)/2 + VIq cos(wot)(cos 5crz+ sniScTx), 



where = ^ a-^^d = 10(^1 + V){A- 

general, ilo is proportional to the slowly varying elec- 
tric field amplitude and the electric dipole moment. 
The time independent part of Hq is diagonalised in 
the energy eigenbasis, — 0(Tf/2, where erf = |5)(.9| ^ 
|e)(e| = sin^a, +cos0a„ {|5>,|e)} = {sin(0/2)|r) + 
cos(6l/2)|0 cos(6>/2 )|r) - sin(6l/2)|0}, 6 = arctan(A/e) 
and (j) = ve^ + A^. We transform Hq to a frame rotat- 
ing at the drive frequency according to 6"*"°"^^*/^, and 
make a rotating wave approximation (RWA), discarding 
terms rotating at frequencies (j) + loq ^ ri,Q, 

HBWA = -{ri<Jt + nal)/2. (1) 

where rj ~ (j) ~ loq, ^ ~ llosinf?, 9 = 6 — 
(5, and ~ denotes the rotating basis. This Hamilto- 
nian is diagonalised in the dressed basis {[— ),[+)} = 
{cos((^/2)|.g) - sin(v5/2)|e),sin(v?/2)|.g) + cos(<^/2)|e)}: 
HimA = -f^Vf/2 where = |-)(-| - |+)(+| = 
siiK^CT^ + cosip = &Yci&n{Q. / rf) £ [0,7r] and il' = 

(ri^ -{-rfY^'^ is the effective Rabi frequency. Note that 
we work in units where H has dimension rad/s. 

We include the influence of phonons, for which the bare 
Hamiltonian is ifph = Sq^q'^q'^q- The electron- phonon 
coupling is generically given by [21] 

i/c-p = 51 *^^q^(q) («q - «-q) ' (2) 
q 

where ^(q) = Ej j'(/ d^r?A*(r)V'i'(i-)e"'''-'")4cj' is the 
Fourier transform of the electron density operator and 
llV'j)} form a discrete basis. Mq = Cq(fi./2/^T^Li;q)^/^ is 
the coupling strength, where is the mass density, V is 
the volume of the lattice and for peizoelectric coupling 
Cq = P, whilst for deformation coupling, Cq = Dq, 
where q = |q|. In the two-level approximation the 
electronic subspace is spanned by the localised states 
{\L),\R)} satisfying (r|L) = V(r - d/2) and {r\R) = 
ip{r + d/2), where d is displacement between dots. Eval- 
uating the integral in g, we discard off-diagonal elements 
such as Mq{L\g{r)\R), since both |Mq| < and {L\R) < 
1 so it is doubly small [1, 3]. The diagonal elements 
evaluate to {L\g{q)\L) = e-*i-'*/2p(q) and (i?|^(q)|i?) = 



2 



giq.d/2p(q-)^ where p(q) = J d^r\iP{r)\^e-''i-r Therefore 
^(q) = p(q)(cos(q.d/2)I — isin{c[.d/2)az), which estab- 
hshes that phonons couple primarily to the position of 
the qubit. For localised states confined to a region of 
length I, the form factor p(q) « 1 for q l/l. The term 
proportional to I just perturbs phonon energies, so 



(b) 



(3) 



where = A/qp(q) sin(q.d/2)/\/y. The Hamiltonian 
for the complete system is i?tot = -f^RWA + ^^o-p + ^^ph- 
We now transform to an interaction picture defined 
with respect to Hhcc = -ffRWA + ^^ph, in which aq(t) = 
aqg-"^"* and af{t) = a^e'^'K Therefore 

Czit) = (cos0cos<y9 — sin0sin(^cos([JoO) '''z 

— (cos 9am(f + sin 9 cos ip cos(woi) ) erf (t) 
~ am9sm{ujot)cryit), 

P^'e-''^'' +Pl,e"^'K (4) 



— cos a sm Lp cr_ , 
- — sin 6* sin (/3 (7^/2. 
Hamiltonian is 



where Pq = cos6'cosiy9crf/2, Pji/ = 
P<^o±iv = Tsin6'(l±cosv3) cr|/2, P^^ 
Then, the interaction picture 
Hi{t) = <Jz{t) Eq5q«e-<>* + flqe-'"-*). 

Master equation. To solve the dynamics of the qubit, 
we develop a master equation for the qubit density ma- 
trix, p. Following [22] we integrate the von Neumann 
equation for the joint density matrix, W ^ of the total 
system, then trace over phonon modes, resulting in 

pi{t) = - f dt' Trph{[i?/(0, [Hi{t'), Wiit')]]}. (5) 

Jto 

We make a Born-Markov approximation, setting t^ = 
— oo and replacing Wi{t') Wi{t), which is valid for 
weak coupling and rapid bath relaxation [22]. When 
and ri', appearing as rotating terms in Hi{t), are much 
larger than the phonon coupling strength, we make a 
second RWA [23, 24], giving the master equation 



PI 



J{uj') {N{lo') + 1) V[P^,]pi + N{u:')V[Pl]pi 



where J(tj) = 271^^ |(jiq|^(5(ijj — Wq) is the spectral den- 
sity, N[uj) = [e.'^I^BT _ j^-j-i -g i-j^g occupation of phonon 
modes, and V\A\p = ApA^ - [A^ Ap + pA^ A)/2. The dif- 
ferent operators P^ are classified according to their effect 
on the system: those proportional to trf produce pure de- 
phasing in the dressed basis, whilst those proportional to 
(7^ induce transitions between the dressed states. 

The second RWA used to derive the master equation 
implies that our analysis is valid only in certain lim- 
its. Firstly, the RWA made above is reasonable when 
u;o,f2' ^ Ji^') so we are considering strong driving and 
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FIG. 1: Energy level diagrams for dressed state raising and 
lowering processes in (a) the dressed basis and (b) the bare 
eigenbasis, for i] < 0, where j— ) ~ |e) and |-|-) ~ \g). Solid ar- 
rows indicate photon absorption, and dotted arrows indicate 
phonon emission. When I77I S> SI' ~ \g\ = udq — (j). 



weak coupling^ where the Rabi frequency is larger than 
the dissipation rate. Weak driving is discussed in [9]. 
This implies that on resonance the DC conductance peak 
is saturated, i.e. the average polarisation of the qubit is 
zero. Secondly, the RWA we have made limits analysis of 
dynamics to times ^ . In what follows we calculate 
steady state properties, so the RWA is not restrictive. 

Setting p(t) = {l+Xd{t)ai+yd{t)(T'l+Zd{t)ai)/2,yie\ds 
the equations of motion for the components of the Bloch 
vector in the dressed basis, id = —{2Tz + {T~+'^+)/2)xd 
(and Xd Vd), and Zd = (r_ -r+) - (F_ +F+)zd, where 



F. = 



( J(0) cos^ 9 cos^ ip(N{{)) + 1 /2) /2 

+ J(wo) sin^ 9 sin^ if{N{uJo) + l/2))/2, 

J(c.o - n') sin^ e{l - cos^f ^^^°~^^^ + ^ 
+ J(c.o + n') sin' 9{1 + cos ^)2 ^K + O0 + ^ 



+J{n') cos' 9 sin' (p{N{n') + (1± l)/2). 



(6) 



Tz contributes to the dephasing rate, but not to transi- 
tion rates between dressed states. The three terms ap- 
pearing in Fzp correspond to interactions with phonons 
of energies luq ± fl' and Q' . For fc^T <^ fl, N and 
the term proportional to J{luq — D,') corresponds to a 
dressed-state raising process, whilst the other two are 
dressed-state lowering processes, shown in Fig. la. The 
steady state of the Bloch vector is Xd = yd = and 
Zd = (r_ — F+)/(r_ + r+), so Zd is determined by the 
balance of the dressed-state raising and lowering rates. 

We now show that driving can induce population in- 
version. For the purpose of discussion we take T = and 
e > 0, and we consider the dynamics below resonance 
77 < 0. For clarity, we describe the most interesting limit, 
1 77 1 3> r2, in which the lower dressed state and the bare ex- 
cited state are approximately the same, |— ) ~ |e). In this 
limit, there is a population inversion wherever dressed- 
state lowering processes dominate, r_ > r+. In this 
limit il' [r^l, sin(p w —R and cos Lp w — l + P^/2, where 
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R = n/n'. Thus cjo - fl' w and ujq + fl' ^ 2wo ~ (j) so 
r+ « J{(j)) si-n^e, 

r_ « J(|77|)cos2 6'i?2 + J(2cJo-(/')sin2 6'i?Vl6• 
Typically, r_ <C r+, since the phonon density of states is 
smah at low energies, so J{\r]\) J{(j)), and 1 away 

from resonance. However, the condition for inversion, 
r_ > F-)-, can be met in certain regimes, namely when 
sin^ 6 is sufficiently small. 

The terms in T± may be understood physically from 
Fig. lb, where the dresscd-state raising and lowering pro- 
cesses at rates T± are shown in the undriven eigcnba- 
sis. In this picture, a dressed-state raising process corre- 
sponds simply to phonon emission accompanied by relax- 
ation of the qubit. The phonon emission rate contributes 
the factor J(0) to F-|., whilst the factor sin^ 9 is the dipole 
matrix element for the phonon-induced transition. In the 
absence of driving this is the only significant process [23] . 

In the presence of driving, dressed-state lowering pro- 
cesses also occur, due to absorption of cither one or two 
photons, accompanied by emission of a phonon of energy 
57' w \ri\ and 2ujq — (j) respectively. Both of these pro- 
cesses contribute to F_, and excite the 2LS. As in the 
the Jaynes-Cummings model of a driven atom [25], away 
from resonance each photon absorption contributes a fac- 
tor R^, giving the factors R^ and i?^ respectively. Each 
photon contribution already implicitly contains a dipole 
transition matrix element. Phonon emission contributes 
factors J(|r?|) or J{2ujo — (j)) respectively. 

Typically the two-photon process is very weak, so in- 
version occurs when J{\j]\)Rl^ cos^ 9 > J {(p) sin^ 9 , i.e. 
when J {\ri\) / J {(f)) > iav? 9 / VP' . Clearly, this inversion 
condition is satisfied for sufficiently small 77 if J is sub- 
quadratic (i.e. J{u}) < auj'^, for some a). These findings 
arc our principle result. Later we show that inversion 
does not depend on the limit \r]\ 3> ^ above, which sim- 
ply leads to a useful physical interpretation. 

The degree of inversion depends on the relative size 
of and J{4>), as well as the temperature. As T 

increases, the phonon modes become thermally occupied, 
so when khT > fl' the thermally-activated, dressed-state 
raising process becomes significant. This is the inverse 
of the second process depicted in Figs, la and lb, where 
relaxation of the qubit from |e) to \g) accompanied by 
stimulated emission of a photon, reducing the inversion. 

We now show that the condition J{ll!) < auP' is sat- 
isfied for piezoelectric phonon coupling (approximately 
constant or Ohmic, depending on phonon dimensional- 
ity), and that this is stronger than deformation coupling. 
If the spectral density were strongly peaked, e.g. due to a 
discrete phonon spectrum, then the two-photon process 
may become significant, analogous to [15]. 

The spectral density J(uj) varies signicantly over the 
frequency ranges of interest and depends on both the di- 
mension of the phonon bath and nature of the coupling 
to acoustic phonons, either deformation or piezoelectric. 





Deformation 


Piezoelectric 


2D 


^(^)^(l-Jo[2.^]) 




3D 


D(;S;)'(l-sinc[2^^]) 


P^(l -sinc[27r^]) 



TABLE I: Spectral density, J{uj), for deformation and piezo- 
electric coupling in isotropic 2D and 3D media, assuming 
p(q) = 1. D and P are the deformation and piezoelectric 
constants, d = |d|, Cs is the speed of sound, /X2 is the areal 
mass density, ^ is the volumetric density, D = 2n'' D^h/ fj,c^d^ , 
P = P^h/2fj,c'^d and Ljp = 2-KCs/d is the angular frequency of 
wavelength-matched phonons. 



We assume the phonon bath is either 2 or 3 dimensional, 
isotropic and dispersionless (wq = Cgq), and the results 
are summarised in Table I. Coupling of GaAs double- 
dots to phonons has been observed in transport measure- 
ments [3], where it was found that spontaneous emission 
rates are dominated by piezoelectric coupling, with the 
effective dimension dependent on the qubit splitting. The 
angular form factor, F , that appears in the 2D spectral 
densities behaves qualitatively like the corresponding fac- 
tor 1 — sine X in the 3D version. Other effects of phonon 
coupling to GaAs dots are discussed in [1, 26] 

To quantify the phonon induced decoherence rates, we 
estimate the characteristic time scales tq = 2mjj'^/DLo'^ 
and Tp = 27ru;p/Pti>o for 3D phonons, from Table I. For 
GaAs, hD = 13.7 cV, hP = 1.45 cV/nm, /i = 5300 



kg/m^ 



5200 m/s [27] and we take d = 500 nm. 



tJo/27r = 24 GHz [6], so hP = 2.3/ieV, ^iD = 33 neV 
and ojq/lop = 2.3. Then tq = 10 ns and rp = 760 ps, 
so piezoelectric coupling dominates, in agreement with 
experiment [3, 5, 6]. 

The variation in the DC response of an electrometer 
is proportional to the time-averaged expectation of the 
double-dot polarisation, z, [23, 24, 28]. From Eq. (4) we 
find z = (cFz) = cos9 cosip Zd- For weak driving, z w e/cf) 
away from resonance, whilst on resonance (77 = 0) z = 0, 
so there is a resonant peak in the electrometer response 
[6, 9]. When |e| » A, the eigenstates are approximately 
localised, e.g. for e » A > 0, we have \g) ~ \l), so 
inversion corresponds closely to localisation in |r). 

Using the rates in Eq. (6), Fig. 2a shows the left dot 
occupation, M = (1 -I- z)/2, versus bias, e, for different 
driving amplitudes, fJo, assuming 3D piezoelectric cou- 
pling. We focus on e > 0, since the opposite regime is the 
same (up to reflections) . There is clear population inver- 
sion, where M crosses 0.5. There are several points to 
note from this figure. Firstly, the population at resonance 
always remains at 0.5. When resonance is approached the 
dressed states have equal time-averaged projections onto 
the localised states, as the qubit performs Rabi oscilla- 
tions through the poles of the Bloch sphere. Hence the 
average polarisation vanishes. Secondly, the peak in the 
conductance shifts to lower energies, as the photon ab- 
sorption probability R^ increases with increasing driving. 
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FIG. 2: (a) Occupation of left dot, M, for increasing Rabi 
frequency from Qq = 0.01 cjo to 0.17ti^o in steps of 0.02 tJo- 
Other parameters are Up = 0.2 uo — A, fcsT = 0.05 tJo and 
S — n/2 corresponding to driving the tunneUing rate, (b) 
Increasing temperature from ksT = to 0.3 a;o in steps of 
0.06 Wo- As above, uip = 0.2 a>o = A whilst flo = 0.16 wo- 



Both of these features arc consistent with recent exper- 
imental observations [14]. Finally, the ripples superim- 
posed on the peak are due to the structure of the spectral 
density, and are spaced by approximately ujp. The promi- 
nence of the ripples follows from neglecting anisotropy in 
the piezoelectric coupling and Cg. In practise, anisotropy 
in Cs (as in GaAs) should suppress the ripples. 

The inversion becomes more pronounced for a 2D 
piezoelectric coupling, since the spectral density for this 
case tends towards a constant for uj > cOp. Transport 
studies of a similar system suggest that the actual phonon 
coupling is intermediate between 2D and 3D [3], so it is 
plausible that experimental results will find a more pro- 
nounced inversion than that shown above. 

In Fig. 2b we show the dependence of M on temper- 
ature. Clearly the resonance is suppressed, and the av- 
erage occupation tends to 50%, as expected. Also, as 
T increases, the peak moves towards smaller e. As dis- 
cussed above, this occurs because when ksT > H.', the 
thermally-activated dressed-state raising process (corre- 
sponding to relaxation from |e) to \g)) becomes signifi- 
cant, which reduces the inversion. Since Q' increases with 
detuning, the inversion is more robust further from the 
resonance condition. This temperature dependence of 
the peak, along with its dependence on driving strength, 
is a signature of the population inversion discussed here. 

Since we are considering a regime of relatively strong 
driving, it is natural to ask whether higher order driving 
terms are significant. Such terms produce harmonic reso- 
nances, and have been seen experimentally [6, 14]. There 
are two origins of harmonic resonances in this kind of sys- 
tem: firstly, at sufficiently high intensity, multi-photon 
processes may become important. Secondly, nonlincarity 
in the coupling of A to the driving voltage, V{t), will pro- 
duce harmonics of loq, which can drive Rabi oscillations. 
Whilst these effects will produce features at cf) = mujo, 
unless f2 ~ Wo [29], they will not significantly affect the 
system near the fundamental frequency. 

In conclusion, we have shown that population inver- 
sion is possible for a driven, dissipative 2LS with sub- 
quadratic spectral density, consistent with recent experi- 
mental observations. The inversion grows as the driving 



increases, and the maximum inversion moves to lower 
energies as both driving amplitude and temperature in- 
crease. This phenomenon is quite generic for physically 
reasonable systems, and does not rely on any specific 
structure in the spectral density. 

We thank Jason Petta, Alex Johnson and Charles Mar- 
cus for previews of experimental data and discussions, 
and Mikhail Lukin and Tim Ralph for useful conversa- 
tions. This work was funded by Fujitsu-CMI and the 
E.U. NANOMAGIQC project (IST-2001-33186). 

Note added: after acceptance at PRL we became 
aware of [30], which also discusses population inversion 
in driven impurity 2LSs. Although the mechanism dis- 
cussed in [30] is similar to that presented here, we have 
determined conditions under which unstructured spectral 
densities produce inversion, and shown that these may be 
satisfied in quantum-dot systems. We thank Mark Dyk- 
man for bringing [30] to our attention. 
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